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^ Abstract: The dynamics of the fluid fields in a large class of causal dissipative 

\ fluid theories is studied. It is shown that the physical fluid states in these theo- 

Q\ • ries must relax (on a time scale that is characteristic of the microscopic particle 

interactions) to ones that are essentially indistinguishable from the simple rel- 
ativistic Navier-Stokes descriptions of these states. Thus, for example, in the 
relaxed form of a physical fluid state the stress energy tensor is in effect indis- 
' tinguishable from a perfect fluid stress tensor plus small dissipative corrections 

qq . proportional to the shear of the fluid velocity, the gradient of the temperature, 

CN ; etc. 

> ' 

oo ; 

^ ■ §1 Introduction 

OO ' 

A simple mathematical model provides an elegant and accurate description of the 
common materials called fluids. The effects of internal dissipation in these materials — 
viscosity and thermal conductivity — are also well modeled by a simple generalization of 
the basic theory called the Navier-Stokes equations. Unfortunately, the most straightfor- 
i i ward approaches to constructing relativistic generalizations of the Navier-Stokes equations 
result in rather pathological theories (Eckart [1], Landau and Lifschitz [2]). These theories 
are non-causal, unstable, and without a well posed initial value formulation (see for ex- 
ample Hiscock and Lindblom [3]). Less straightforward approaches have succeeded more 
recently in producing a class of causal dissipative fluid theories (e.g., Israel and Stewart 
[4], Carter [5], Liu, Miiller, and Ruggeri [6], Geroch and Lindblom [7], etc.). These the- 
ories have eliminated the pathologies of the straightforward relativistic generalizations of 
the Navier-Stokes equations, but they do so at the expense of increasing significantly the 
number of dynamical fields needed to describe the fluid. Unfortunately the additional 
dynamical degrees of freedom associated with these extra fields have never been directly 
observed in real fluids. This is probably why these new theories have not found widespread 
acceptance. 

In this paper the dynamics associated with these additional fluid fields are analyzed 
in a very large class of causal dissipative fluid theories. It is shown that the physical fluid 
states relax (on a time scale characteristic of the inter-particle interactions) to ones that 
are also well described by the simple relativistic Navier-Stokes theory. For example, the 
stress-energy tensor in such a relaxed fluid state is well described by the usual perfect 
fluid stress-energy tensor plus the Navier-Stokes expressions for the dissipative corrections 
involving the shear of the fluid velocity, the gradient of the temperature, etc. This result 
suggests that meaningful differences between the causal theories and the non-causal Navier- 
Stokes theory can not be observed. The complicated dynamical structure of the causal 
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theories is necessary to insure that the fluid evolves in a causal and stable way. But this 
rich dynamical structure is unobservable, since the physical states of a fluid always evolve 
in a way that is also well described by the Navier- Stokes expressions for the stress-energy 
tensor, etc. The arguments which lead to these conclusions are extremely general: they are 
based on a fully non-linear analysis of the equations and do not assume that the fluid state 
is close to equilibrium. This analysis generalizes significantly the studies of the analogous 
relaxation effect in the solutions of the hyperbolic heat equation (see Nagy, Ortiz, and 
Reula [8]), and the studies of the relationship between the relativistic Navier-Stokes and 
the causal fluid theories in the near equilibrium fluid states (see Hiscock and Lindblom 
[9])- 

Let us begin by recalling the theory of a perfect fluid: the mathematical description 
of a fluid having negligible internal dissipation. The state of such a fluid is determined by 
three fields on spacetime: a future directed unit timelike vector field, u a , and two scalar 
fields n and p. These fields are assumed to be solutions of the differential equations 

V m N m = 0, (1) 
V m T ma = 0, (2) 

where N a and T ab are given in terms of the fluid fields by 

N a = nu a , (3) 

T ab = (p + p)u a u b + pg ab . (4) 

Here p is a smooth function of n and p (the equation of state), that is fixed once and 
for all for a given type of fluid. The conserved vector N a is the particle current of the 
fluid, and thus u a may be identified as the four-velocity and n as the number density as 
measured by an observer co- moving with the fluid. The conserved tensor T ab is the stress 
energy of fluid. Thus from eq. (4), p is identified as the mass-energy density and p as 
the pressure of the fluid, both as measured by a co-moving observer. These quantities 
are all directly observable because the particle current N a and the stress energy T ab are 
themselves directly observable. 

The theory of a perfect fluid, eqs. (l)-(4), has a number of attractive mathematical 
properties. One of the most important of these is that eqs. (l)-(2) form a symmetric- 
hyperbolic and causal system when suitable restrictions are placed on the equation of state. 
Let £ a = (n, p, u a ) denote the dynamical fluid fields. Then eqs. (l)-(2) are equivalent to 

M m a(3 V m ^ = 0, (5) 

where 

d N m Qj> am 

The quantities P a and P aa (functions of £ a and the spacetime metric g a ^) may be chosen 
so that these equations are symmetric in the sense that M m a p = M m p a (see Ruggeri and 
Strumia [10], Geroch and Lindblom [7], and §111 below). When suitable restrictions are 
placed on the equation of state these equations are also hyperbolic and causal, because 
X m = M m a pZ a zP is past directed timelike for every choice of Z a ^ in these theories. 

Next, let us turn to the main subject of this paper: theories of dissipative fluids. 
Since there is as yet no universally accepted theory for such fluids, a rather broad class 
of theories has been included in this discussion. The state of the fluid in these theories is 
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determined by two sets of fields, £ a and (p , each representing some collection of tensor 
fields (possibly subject to certain algebraic constraints) on spacetime. The £ a are to 

represent, as in the perfect fluid case, the dynamical fluid fields. The cp are to represent 
additional 'dissipation' fields that are needed to complete an acceptable causal fluid theory. 
It seems reasonable to restrict the dimension of the combined £ a and p A spaces to be 
equal to the number of independent observable fields in the theory. In the case of a 
simple dissipative fluid — the case that is of primary interest here — the particle current iV a 
and the stress-energy tensor T ab are the independent observable fields. Hence the most 
appropriate choice for the dimension of these combined spaces is fourteen in this case. For 
now, however, neither the structures nor the dimensions of these spaces will be restricted. 
The fields, £ a and (p A , are assumed to be solutions of the system of equations 

M m ap V m ^ + M m aA V m p A = 0, (7) 

M m AB V<p B + M m aA VU a = -lABV B - (8) 

The quantities M m a p, M m aA , M m AB , and I AB are assumed to be smooth functions, 

fixed once and for all for a given theory, of the fields £ a , cp A , and the spacetime metric g a i,. 
Thus eqs. (7)-(8) form a first-order system of partial differential equations for the fluid 

fields £ a and cp A . Three conditions are now imposed on this system of equations. These 
conditions are very general and should apply to essentially any theory of fluids (including 
those describing superfluids, mixtures of different kinds of fluids, etc.). 

Condition i) The first condition is on the M's that appear on the left sides of eqs. (7)- 
(8). Assume that the M's are symmetric, M m a(} = and M m AB = M m ^ AB y, and 

assume that every vector A m given by 

A m = M m a/3Z a z + 2 M m aA Z a Z A + M m AB Z A Z B , (9) 

for some (Z a , Z A ) ^ is past-directed timelike. This is just the condition needed to insure 
that the system (7)-(8) is symmetric, hyperbolic, and causal (see for example Geroch and 
Lindblom [11] or Miiller and Ruggeri [12]). 

Condition ii) The second condition involves the tensor I AB that appears on the right side 

of eq. (8). Assume that I AB Z A Z B > for every Z A ^ 0. 2 This condition is adopted to 
insure, as will be seen more clearly below, that this fluid theory is strictly dissipative. It is 
precisely analogous to requiring that the viscosity coefficients and the thermal conductivity 
not vanish in the Navier-Stokes equation. 

Condition Hi) The third condition concerns the conservation laws. Assume that there 

exist specific smooth functions A^ a and T ab of the fields £ a , cp A , and g a ^ such that eq. (7) 
implies the conservation laws, eqs. (1) and (2). This condition merely insures that the 
theory possesses a conserved stress energy tensor and particle current. 



This restriction is not required in the analysis presented here however. If the number of 
fluid fields were taken to be larger than the number of observables then some of the fluid 
fields would not be observable. The results derived here would still apply, but some of 
them would change character from experimentally testable predictions to mathematical 
identities. 
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The main result of this paper is derived in §11. It is shown that physical states of the 
fluid relax — on a time scale r that is characteristic of the inter-particle interactions — to 
ones in which the dissipation field is determined in effect by the dynamical field £ a and its 
derivative. In particular, a bound is derived for the quantity Acp^, defined by 



(10) 



</ = - \{i- l ) AB M m aB v m e] + A</, 

in the physical fluid states of any fluid theory which satisfies Conditions i)-iii). This bound 

on Aip A is smaller by the factor 6 5 / 2 (rf /L) 2 than it is expected to be. The constant v is a 
characteristic sound speed, L is a macroscopic length scale that characterizes the particular 
state of the fluid, and the constant b is a dimensionless bound on M m ^g, M m a A, [I~^) AB 
and their derivatives (which will be defined precisely in §11). The constant b is expected 

to be of order unity for 'reasonable' fluid theories. Thus the factor 6 5 / 2 (rf/L) 2 should be 
extremely small for real fluids. For example in water (tv/L) 2 ps 10~ 12 for fluid states with 
L 0.1cm. Since the dissipation field in a relaxed fluid state is determined in effect by 
the dynamical field £ a and its derivative, then so are all other functions of the fluid fields. 
In particular, the particle current N a and stress energy tensor T a ^ are given by 



N a = 



tab 



3N a / t\AB 

Na -^i ( J ) Mm *Bv m e 



tp c =0 



+ AN a , 



T ab _ dT 



ab 



dip" 



AB 



M m aB v m e 



+ AT 



ab 



(11) 



(12) 



ip c =0 



It is shown that the quantities AN a and AT a ^ are also smaller than their expected values 
by the factor b^(rv/L) 2 . These results apply to any dissipative fluid theory that satisfies 
Conditions i)-iii) above, and to any physical fluid state (i.e., as defined more precisely 
below, a state in which the spatial and temporal variations of the fluid fields are larger 
than the microscopic scales). This result explains why the independent dynamics of the 

dissipation field cp A is never observed: its value is determined in effect by the dynami- 
cal field £ a and its derivative, via eq. (10), on any time scale over which a macroscopic 
measurement of the system can be made. Although measurements could in principle be 
carried out on fluid systems over very short time and distance scales, it is not required or 
even expected that such measurements will be modeled in detail by any macroscopic fluid 
theory. 

The results of §11 show that a dissipative fluid quickly relaxes to a state in which 
the particle current and stress-energy tensor are determined (in effect) by the dynamical 
fluid field £ a and its derivative V TO £ a . Such relationships are quite familiar to us; for these 
are precisely the forms that the expressions for these quantities take in the Navier-Stokes 
theory. Recall that in the relativistic Navier-Stokes theory (as formulated by Eckart [1]) 
the particle current and stress-energy tensor are given in terms of the fields £ a = (n, p, u a ) 
by 



2 The somewhat more general function — /^(^ a , <p B ', g a ^) could have been adopted for the 
right side of eq. (8) if it satisfied a few additional constraints. This more general form is 
equivalent to that given in eq. (8) if and only if 7^ satisfies the following three conditions: 



a) I a = when <p u = 0, b) I Ay > when ip B ^ 0, and c) dlj^/dip is not degenerate 
when evaluated at cp C = 0. 
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(13) 



(p + p)u a u + pg' 



ab + T ab + Tq ab + 2u (a q b)^ 



(14) 



where 




— Or, r ,am rl bc ^ rl ab •cm v7 

- Zr]i q q - -q q V (jn u c) 



(15) 



m 



(16) 




« (q am V m T + Tu m V m u a ) . 



(17) 



The quantities 771 , 772, and k (positive functions of n and p) are the viscosities and thermal 
conductivity respectively; and, the quantity T (a function of n and p) is the thermodynamic 
temperature which satisfies the first law of thermodynamics, 



In this theory the conservation laws, eqs. (l)-(2), and eqs. (15)-(17) are the differen- 
tial equations that determine the field £ a . If the dissipation fields are defined as, <p A = 
(r, q a , T a ^), then this theory is of the same general form as those being considered here. 
The conservation laws are precisely in the form of eq. (7), while eqs. (15)-(17) have the 
form of eq. (8). This theory fails to be an acceptable theory because M m j^q = and 
thus it fails to satisfy Condition i). Note that for this relativistic Navier-Stokes theory, the 

quantities Aip A , AN a , and AT ab as defined in eqs. (10)-(12) vanish identically. 

The vanishing (effectively) of AN a and AT ab for the general dissipative fluid theories 
considered here implies that the particle current and stress-energy tensor depend (in effect) 
only on £ a and its derivative V m £ a in any physical fluid state. In the relativistic Navier- 
Stokes theory, however, only certain components of V TO £ a appear in these expressions. 
For example, in the Navier-Stokes theory V m T appears in these expressions but not the 
gradient of any other thermodynamic scalar. It is natural to ask then, what class of 
fluid theories have the property that their fluid states always relax to ones in which the 
particle current and stress-energy tensor are in effect indistinguishable from those of the 
relativistic Navier-Stokes theory? Or in particular, in which theories do N a and T a ^ depend 
on £ a and V TO £ a in precisely the same way as in the Navier-Stokes theory? The following 
two additional conditions are necessary and sufficient to guarantee that a theory will be 
indistinguishable from Navier-Stokes in this way: 

Condition iv) The fourth condition concerns the space of dynamical fields £ a . Assume 
that eq. (7) is precisely equivalent to the conservation laws, eqs. (l)-(2). This implies that 
the space of the £ a consists of one vector and one scalar field which may, without loss 
of generality (as shown in §111), be taken to be £ a = (n, p, u a ), where nu a = N a (with 

u a u a = —1) and p = u a iii ) T ab . 

Condition v) The fifth condition concerns the tensor M m a A that appears in eqs. (7)-(8). 
Assume that M m ' m £, a depends on V TO n and V m p only in the combination V m T = 
{dT/dn) p V m n + {dT/dp) n V m p in the (p A = states of the fluid. This condition is required 
to insure that heat flow is generated by the gradient of the thermodynamic temperature 



dp = nTds + 



p + p 



dn. 



(18) 



n 
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T and not the gradient of some other thermodynamic scalar. This condition is equivalent 
to the requirement that the equilibrium states of the fluid be 'isothermal.' 

The theories that satisfy these two additional conditions are the natural causal generaliza- 
tions of the Navier-Stokes theory: the causal theories of a simple dissipative fluid. 

In §111 the expressions for the relaxed forms of the particle current and stress-energy 
tensor are evaluated for the theories of simple dissipative fluids, i.e., those satisfying Con- 
ditions i)-v). With the convenient choice of dynamical fields, £ a = (n, p, u a ), eqs. (1.11) 
(1.12) reduce to 



N a = nu c 



(19) 



T ab =(p + p)u a u b +pg ab + 2 m 



q am q bc 



l q ab q cm 



V {m u c) + m V m u m q ab 



2 KU ( a 



q b ^ m V m T + Tu\ m \v m u b ^ 



+ AT 



ab 



(20) 



for suitably chosen functions (of n and p) p, rji, 772 , and k. Since AT ab is extremely small 
in the physical fluid states of these theories, this shows that the particle current and stress 
energy tensor are (in effect) indistinguishable from those of the relativistic Navier-Stokes 
theory. 3 



§11 The Relaxation Effect 

The key result in this paper is that the physical states of the fluid relax to ones in 
which the dissipation field cp A is determined (in effect) by the dynamical fluid field £ a 
and its derivative V m £ a . That some form of relaxation should occur in the solutions of 
eqs. (7)-(8) can be seen fairly easily. Consider the quantity I ab<£ B ' + M m 'm£ Q • If this 
quantity does not vanish at some point, then the first term in eq. (8) causes <p A to evolve 
in the direction that tends to make it vanish. The rate at which this evolution occurs is 
determined by the time scale that is encoded in the tensor I A B- F° r fluids this time scale 
will be determined by the viscosity and thermal conductivity coefficients contained in Ij^Bi 
and therefore will be characteristic of the inter-particle interaction times for the fluid. The 
demonstration that the quantity Acp^ defined in eq. (10) is small will be done in two steps. 

First, it is shown that a related quantity cp A + A A , defined below, is small using a fairly 
simple and straightforward argument. Second, a slightly more elaborate argument shows 
that the quantities A(p A , AN a , and AT ab of eqs. (10)-(12) are also small. 
Begin by obtaining the following equation for ip A + A A from eq. (8): 



m 



= -21 AB (V + A A ) (<p B + A B ) + (V + A A ) A A , 



(21) 



Note that AN m and UaU^AT ^ vanish identically as a consequence of the particular choice 
of £ a made here. Had a different choice been made, such as the one traditionally used in 
the Landau-Lifschitz theory [2], then other components of these quantities would have 
vanished identically instead. 
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where K A and A A are denned by 

h A = (V 1 ) AB [M m aB V m e ~ \^ m M m BC ^ , (22) 

A A = A B V m M m AB + 2M m AB V m A B . (23) 

Next consider S(0), a bounded open subset of some Cauchy surface. Use the timelike 
vector field whose divergence appears on the left side of eq. (21) to define a map between 
the points on successive Cauchy surfaces. Let S(k) denote the image of S(0) under this 
map into the Cauchy surface labeled by the time function k. Choose this time function k 
so that it satisfies 

I AB Z A Z B > -Z A Z B M m AB V m n (24) 

(for every Z A ), in the spacetime region (0, k ) x S(k ). Next, define the following C? norm 
of^ A + A A , ' 

« 2 («) = J ) Gab (/ + A A ) + A S ) dV, (25) 

where Gab = n mM m AB , and n m is the future directed unit vector proportional to V m K. 
The evolution of this norm is determined by integrating eq. (21) over the spacetime region 
consisting of points in S(k) that lie between two nearby k = constant slices. The integral 
along the timelike boundary of this region vanishes because of the choice of S(k). The 
integral of the terms on the right in eq. (21) may be transformed using eq. (24) for the 
first term and the Schwartz inequality for the second. Taking the limit as the difference 
between k on these two slices goes to zero, the following differential inequality is obtained 
for a, 

|<-« + 5lWI, (26) 

where 



\\A\\(k) 



G AB A A A B dV 



S{k) -Vm 



(27) 



and G AB denotes the inverse of G A B- This ordinary differential inequality, eq. (26), can 
be integrated to obtain the following bound on a, 

a(K ) < a(0)e- K ° + J e -( K °-^\\A\\(K)dK. (28) 

2 JO 

To proceed further a bound must be obtained for the quantity ||«4|| that appears 
in eq. (28). To this end a norm is introduced on tensors: The positive definite G a B = 

n m M m a Q and its inverse G a ^ are used for indices associated with the dynamical field, 

and the positive definite Gab = n mM m AB and its inverse G AB are used for indices 

associated with the dissipation field (p A . For spacetime indices the positive definite metric 

Gab = n aJi>b + v ~ 2 (9ab + n aTib) and its inverse G ab = n a n b + v 2 (g ab + n a n b ) are used. 
The constant v, with < v < 1, is chosen to be an upper bond on the speed (relative 
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to n a ) of signal propagation, i.e., a number such that (n a A a ) 2 > v 2 (g a \) + n a ni ) )X a X b 
for every A a given in eq. (9). As examples of this norm, the integrand in eq. (25) can 
be written as \^ A + A A \ 2 = Q AB (^ A + A A )(ip B + A B ), while \A A \ 2 = Q AB A A A B and 



lV m «;| z = Q ah W a ^b K = -V a KV a K. Note that \M m AB \ < Ad where d is the dimension of 



the space of dissipation fields. Since d will be some relatively small integer, say d 
norm of the M's will be of order unity in these fluid theories. 



9, the 



In the fluid theories considered here the quantities M m AB , M m aA and I AB are 
assumed to be smooth functions of the fluid fields. Therefore, these quantities and their 
derivatives with respect to the fluid fields are bounded. It is convenient to quantify these 
bounds in terms of three constants 6, r and Consider fluid fields £ a and cp A that are 
bounded by the constant (: 



</]<c, n<c 



(29) 



Next define the dimensionless constant b to be a bound on the M's and their derivatives. 
In particular assume that 



|M 



d 2 M r 



m 



AB\ 



AB 



dM m AB 


b 


dM m AB 


d f.a 
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d 2 M r 



AB 



d£ a d<p c 



d 2 M 7 



AB 



d(p c dcp 



D 



(30) 



and 



M rn aA \<b, 



9M m aA 


b 


dM m aA 






dip B 



< 



c 



(31) 



Finally, the constant r is defined as a bound on (I l ) AB and its derivatives 



\(i-y B \<bT, 



UAB 



br 

£ 7' 



-UAB 



dip 



C 



(32) 



for b given above. The constant r that appears in these bounds is the characteristic time 
scale on which the dissipative term I AB influences the evolution of the fluid in eq. (8). 
This constant also fixes the relationship between physical time and the time function k 
because of eq. (24). The time function k can be chosen so that 



|V TO tt| > 



1 



(33) 



This k in effect measures time in units of r. 

To proceed further bounds must now be placed on the spatial derivatives of the fluid 
fields. Assume that there exists a constant L such that 



|V m £ a | < ^, \V m ^ A \ < ^, |V m V n £ a | < ^, |V m V n </| < (34) 

These inequalities restrict the solutions to the fluid equations 4 to those which do not vary 
appreciably on length scales shorter than L and on time scales shorter than L/v. These 



It is expected that large numbers of solutions to the fluid equations exist which satisfy 
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inequalities select, therefore, the set of solutions that represent real physical fluid states. 
Fluid states in which rapid variations of the fluid fields occur on time and length scales 
smaller than the microscopic particle interaction scales probably can not be adequately 
modeled by any macroscopic fluid theory. Thus, solutions to the fluid equations having 
these properties are not considered physical. The inequalities in eq. (34) therefore select 
out the physical solutions of the fluid equations when L is larger than the microscopic 
interaction length scale. For these solutions the quantity ||^4|| can be bounded by using 
eqs. (30)-(34) in eq. (27): 

||^||(/,)<266 3 C(^)V 1 /2( K ), (35) 
where V(k) is the volume of the region S(k): 

V(k) = [ dV. (36) 
JS(k) 

Including this bound into eq. (28) the following bound is then obtained for a: 

2 C 

a(«o) < a(0)e- K ° + 136 3 C (^) J e'^-^V^i^dK. (37) 

This bound consists of two pieces. The first is simply the initial value of a multiplied by 
e -K This term falls exponentially to zero on the characteristic time scale r. The second 
term is a constant multiplied by the time average of the spatial volume over which the 
norm a is defined. This second term, the asymptotic bound on the norm a, is smaller 
than the a priori expectation of its value, ((bV)^^, by the factor 6^/ 2 (rf /L) 2 . This factor 
will be extremely small, being proportional to the square of the ratio of the characteristic 
dissipation time scale r to the characteristic dynamical time scale L/v, as long as the 
constant b is of order unity and the constant v is comparable to the sound speed in the 
material. The constant b is a measure of the M's and I and their derivatives with respect 
to the fluid fields. This constant will be of order unity unless these quantities depend 
on the fluid fields in a very perverse way (e.g. if the M's dependence on the fields were 
highly oscillatory) . The constant v will be comparable to the sound speed of the material 
as long the foliation of Cauchy surfaces is chosen so that the fluid motion is not highly 
supersonic, and as long as the characteristic speeds associated with the dissipation fields 
are comparable to the usual sound speed. In this case the bound on a derived in eq. (37) 

implies that the dissipation field cp A relaxes in the physical states of these fluid theories in 
such a way that the quantity (p A + A A becomes extremely small. 

To complete the argument that the quantities Atp A , AiV a , and AT ab of eqs. (10)- 

(12) are small, additional C? and £ 4 bounds are needed on the dissipation field (p A . To 
obtain these bounds the following identities are derived from eq. (8), 



m 



M 



"WV] = -Vab<p\ B + <P A B A , (38) 



m 



M m AB g CD <p A <p B <p c <p D 



-2I AB g C D<p\ B tp C tp D + ip A ip B ip C C ABC , (39) 



these conditions. In particular, it is expected that initial data satisfying these conditions 
on a Cauchy surface will evolve (for some macroscopic time) as a solution that satisfies 
these conditions everywhere in the development of these data. There do not exist theorems 
at present, however, which prove the existence of solutions having these properties. 
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where B A and C AB q are defined by 



B A = V B V m M m AB - 2M m aA V m e, 

C ABC = GAB^C + 2M m AB g CD V m ip D + M m AB tp D V m GcD- 
Next an integral norm, analogous to a above, is defined for the field (p A : 



(3 2 {k) 



j. 


~K\ 


V A + A A 


2 9 




2 „ 
+ K% 




4" 


IS{k) 

















dV, 



(40) 
(41) 

(42) 



where K\ > 0, K 2 > 0, and K% > are constants whose values will be specified later. The 
time evolution of j3 is determined in analogy with eq. (26) by integrating K\ multiplied 
by eq. (21), plus K% multiplied by eq. (38), plus multiplied by eq. (39), over the 
spacetime region consisting of points in S(k) that lie between two nearby k = constant 

slices. The sequence of spatial sections S(k) is chosen in this case so that the timelike 
boundary integral vanishes identically here as well. The integrations on the right side of 
this equation may be simplified, again in analogy with eq. (26), by using eq. (24) and the 
Schwartz inequality. The result is the following differential inequality on (3, 



^<-p +1(^1^11 + K 2 \\B\\) + ^K 3 



\Cabc\' 



D 



dV 



S(k) 



|V m /c|' 



1/2 



(43) 



where 



is given by eq. (27) and ||£>|| is 
OT(«) = 



\B A \ 2 dV 



S{k) |V m K| 



1/2 



(44) 



The quantities that appear on the right side of eq. (43) can be bounded if one additional 
restriction is made on the physical solutions of the fluid theory. Assume that the extrinsic 
curvature and acceleration of the k = constant surfaces are bounded by 



|V a n & | < — . 



(45) 



Using the bounds given in eqs. (29)-(33), (34), and (45), the following bounds can be 
obtained for the quantities that appear on the right side of eq. (43): 



|£||<4&Cp 1/2 (^ 



S(k) 



\ C ABC\ 2 




2 

dV 






2 



n l/2 



< 96 2 C— A 
~ L K2 



(46) 



(47) 



Combining these bounds with eq. (35), the differential inequality for (3 can be simplified 
to the following 
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d/3 
d,K 



< - 



1 - h\™^ 
2 L K2 



P + 



l^(^)K l + 2K 2 }bc(^)v^\n). (48) 



The constants K\, K2 and K3 are now chosen to be 

1 



136 2! 



K 2 = -(- 



K 3 = 



186 2 C' 



With these choices eq. (48) becomes 

dp p 



(49) 



(50) 



Integrating this inequality, the desired bound on the norm j3 is obtained: 

PM < (3(0)e' Ko/2 + 46C (y) 2 < >> 
where < V 1 / 2 > denotes the time average of the spatial volume, 

<V 1 / 2 >=- fe-M)^^ 
2 JO 

This bound on (3 implies an C 2 bound on <p A + A^, and simultaneously C 2 and £ 4 bounds 
on 99^. The asymptotic values of these bounds are given by 



(51) 



(52) 



<p A + A M \'dV 



A\2. 



S(n) 



1/2 



>, 



S(k) 



\<P 



A\2 



dV 

1/2 



1/2 



<86C 



TV 



< 



>, 



S(k) 



dV 



<72b\ 2 (^) 2 <^ 1/2 >- 



(53) 



(54) 



(55) 



These bounds on cp A are smaller by the factor rv/L than their a priori expected values. 
The bound on cp^ + is even smaller, however, being reduced from its a priori expected 
value by the factor (tv/L) 2 . Thus (p A + A^ becomes small not simply because tp A and 
A^ get small individually. Rather, this quantity becomes small because cp A approaches 
A^ asymptotically. Note that the region S(k) over which these norms are computed may 
be chosen arbitrarily on any particular slice. 5 Also note that < V 1 / 2 (k) >^ V 1 / 2 (k) if 

V(k) ^> (tv ) 3 . The time average used here is exponentially weighted and hence only those 
slices within about one microscopic interaction time r of k contribute significantly. 



5 The regions S(k) on the other slices in the foliation are then fixed, however, in order to 
eliminate the spatial boundary terms from the integration. 
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The main results of this section are bounds on the quantities A(p A , AiV°, and AT ab 
defined in eqs. (10)-(12). These bounds are obtained beginning with the quantity Atp A , 



A(p A = </ + {r l ) AB M m aB v m e 



tp c =G 



(56) 



This quantity can be re-written as the sum of cp A + A A , a quantity whose bound was 
established above, plus e A : 



A(p A = </ + A A + e A , 



(57) 



where e A may be written (using the standard expression for the remainder in a Taylor 
expansion) as 



A 1 ( T -1\AB,(3y7 njm 



d 



o \d(p c 



(I- 1 ) AB M m aB \ (^,\p D )\d\. 



(58) 



Now, using eqs. (30)-(32) and (34) it is straightforward to obtain the following bound on 



e A 



\e A \ < 36 2 (^) |^|. 



(59) 



Using the triangle inequality for C? norms, the norm of A(p A can be expressed as the sum 
of the norms for cp A + A A , from (53), and the norm of e A , using (54) and (59): 



S(k) 



\A<p A \ 2 dV 



1/2 



< 386 



< v 1 ' 2 > . 



(60) 



Thus, the norm of Atp A is smaller than its a priori expected value by the factor o^I 2 (tv / L) 2 . 
Turn next to the quantity AT ab , 



^rpab rpab 



rpab _ 



dT 



ab 



{i- l ) AB M m aB v m e 



(61) 



This quantity may be re-written (using eq. [56] and again the standard expression for the 
remainder in a Taylor expansion) as 



AT ab = Aip A 



dT ab 
dip A 



<p c =0 



+V ^ B M {l _ x) ^_ {eAtpC))dy 



(62) 



The norm of this quantity can easily be bounded by 



|AT a& |<-|A<^| + ^|^ 

if the field derivatives of T ab satisfy the following bounds 
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(63) 



dT ab 



dip 



A 



e 



2rpab 



d z T 



(64) 



The constant e is a characteristic internal energy density. Using the expressions for the 



bound on Aip from eq. (60) and the bound on \<p | from eq. (55), the following bound 
on AT ab is obtained, 



\AT ab \ 2 dV 



S(k) 



1/2 



< H26 3 r 



TV\< 

V~l) 



e < 



> 



(65) 



This equation provides a bound on AT ab that is smaller than its a priori expected value, 
e < V 1 ! 2 >, by the factor b 3 (rv/L) 2 . 

An exactly analogous bound can be obtained for AN a if the field derivatives of N a 
are bounded by 



dN a 



dip 



A 



v 



2 ATCl 



d z N 



d(p A d(p 



B 



v 



(66) 



where v is a characteristic number density. The bound on AN a is obtained in precisely 
the same way as the bound on AT ab , with the result 



JS{k) 



\AN a \ 2 dV 



1/2 



<1126 3 (y)^<^ 1/2 > 



(67) 



Thus the bound on AJV a is also smaller than its a priori expected value by the factor 
b 3 (rv/L) 2 . 



§111 Simple Dissipative Fluids 

In this section the relaxed expressions for the particle current and stress energy 
tensor, eqs. (11)— (12), are evaluated for the theories of a simple dissipative fluid. Condi- 
tion iv) guarantees that eq. (7) is equivalent to the conservation laws in this case. This 
implies that the space of the £ a must consist of one vector and one scalar field. The 
form of eqs. (7)-(8) is unchanged if the fluid fields are transformed in the following way: 
£a _ ^(^ ipB^ anc [ (pA _ ^A^By r-Q^ c hoice £ a = (n, p, u a ) and (p A = cp A , where 

nu a = N a (^,<p B ,g bc ), (68) 



p = u a u b T a \^,<p»,g cd ), (69) 

is a transformation of this form. Thus £ a may be chosen to be £ a = (n, p, u a ), without loss 

of generality. With this choice dN a /dcp^ = and u a u b dT ab / d(p A = 0. Evaluating eq. (11) 
for this case we obtain N a = nu a + AN a , hence eq. (19). The quantity AN a vanishes 
identically as a consequence of the choice of £ a used here. Condition iv) also implies that 
the tensors M m a p and M m a ^ °f ec L- (?) must be given by 
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8N m Qrpam 
+ Paw 



M 



m 



dT 



am 



aA 



dp' 



(70) 



(71) 



where P a and P aa are suitably chosen functions of £ a , ip^ and g^. Note that the term 
proportional to P a is missing from eq. (71) because dN a /dip^ = for our choice of £ a . 
Using this expression for M m aa , the general expression for T ab in eq. (12) reduces to 



T ab _ dT 



dip 



A 



AB dT 



cm 



dip 



D 



'Pac^m£, a 



+ AT 



ah 



(72) 



ip c =Q 



Condition v) places restrictions on the allowed forms of M m a j\V m ^ a in the fluid 
states where ip^ = 0. From eq. (71) it follows that this quantity is determined by P a a- 
In the ip^ = fluid states the tensor M m a g is identical to the tensor that governs the 
evolution of a perfect fluid via eqs. (5)-(6). The most general P's that make M m a p 
symmetric in this case are given by 



Pade = -Qi 



dp\' 
dn J 



+ Q2 



p + p 



n 



dn — Qi 



dp\ f dp^ _ q 2 P + P 



dp J \dnj 



n 



dp, (73) 



Paadi a = U a Ql 



dp 
d~p 



n 



dp\ q 2 p + p 
dn J p n 



dn + u a Qi 



Qabdu 



dp 
dp 



+ Q2 



n 



dp 



(74) 



where Q\ and Q2 are arbitrary functions of n and p, and g a 5 = g a b + UaU^ (see Geroch 
and Lindblom [7]). The hyperbolicity and causality conditions for M m a p in this case are 
simply, Q\ > and Q2 > 0, and the equation of state must satisfy 



< 



< 1, 



(75) 



dp 

PP; s 

with n > and p + p > 0. Thus the tensors M m a n and M m a ^ are determined completely 

(up to the arbitrary overall factor Q\) in these fluid states by the function Q^- 

Condition v) fixes Q2 by demanding that M m ay (V m ^ a and hence -P a ( a V m )£ a de- 
pend on V TO n and V TO p only in the combination V TO T = (dT / dn) p V m n + (dT/dp)nV m p. 
The unique Q2 which insures this is 



Q2 



1 



dp 
d~p 



&T\ 



dp 
dT 



n 2 T 2 



dp 
dp 



dp 
~d~s 



dp\ 
ITs) 



e 



(76) 



where T and s are the temperature and entropy that satisfy the first law of thermody- 
namics, eq. (18), and 6 = (p + p)/nT — s. The second equality in eq. (76) shows that the 
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condition Q2 > 0, needed to insure hyperbolicity of the equations, is equivalent to a well 
known condition for thermodynamic stability (see Hiscock and Lindblom [13]). With this 
choice of Q2 the quantity -P a ( a V m ^ a reduces to 



P + P 
T 



dp 
dp 



u {aSm) T 



TV 



(to u a) 



(77) 



The tensor dT ah / dip A (I- l ) AB dT cd / dip B that appears in eq. (72) depends only on £ a and 
g ab . The most general such tensor (having the appropriate symmetries, etc.) depending 
only on £ a and g ab is given by 



dT 



ab 



dp 



(- 1 ) 



AB dT cd 



dy B Qi(p + p) \dp 



q a{c q d)b _ \ q ab q cd 
3 



+ m q ab q cd - 2KTu^ a q b ^ c u d ^ 



(78) 



The arbitrary functions 771, 772, and k (of n and p) that appear in eq. (78) must be positive 
as a consequence of the positivity of IaBi horn Condition ii), and the positivity of Q\ 
and (dp/dp) s , from the hyperbolicity and causality of the equations. 6 Combining this 
expression, eq. (78), with eq. (77) in eq. (72) results in the desired form, eq. (20). Thus, 
the relaxed form of the stress energy tensor is indistinguishable in these general causal 
theories from that of the relativistic Navier-Stokes theory. 



§IV Concluding Remarks 

The argument presented here demonstrates that a relaxation effect takes place in 
virtually every causal theory of dissipative fluids. In the relaxed fluid states the stress 
energy tensor and particle current are well described by expressions that depend only on a 
subset of the fluid fields (referred to here as dynamical fluid fields) and their derivatives. For 
those theories that represent simple dissipative fluids, these expressions are identical to the 
ones given by the relativistic Navier-Stokes theory. This implies that any measurement of 
the stress-energy tensor or particle current in these theories (made on any time and length 
scale that exceeds the microscopic particle interaction scales) will give results that are in 
effect indistinguishable from those of the Navier-Stokes theory. Of course the Navier-Stokes 
theory is not really a proper physical theory at all since it is non-causal, unstable, etc. It is 
incapable of predicting the future evolution of initial fluid states. The argument presented 
here shows, nevertheless, that the evolution of any physical fluid state according to any 
causal theory results in stress-energy tensors and particle currents that are experimentally 
indistinguishable from the Navier-Stokes expressions for these quantities. Further, this 
argument shows that the independent dynamics associated with the dissipation fields of 
the fluid (i.e., those additional fluid fields that are added to the theory to make it causal) 

The only requirement on the dissipation fields needed to obtain eq. (78) is that the space 

of the <p A be large enough to insure that none of the coefficients 771, 772, or k vanishes 
identically. This requires in particular that this space be at least as large as the nine- 
dimensional space of symmetric trace-free tensors. This is precisely the dimension that is 
appropriate for a theory in which the particle current N a and stress energy tensor T ab are 
the only independent observable fields. 
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is not directly observable in the physical fluid states. On a time scale that is characteristic 
of the inter-particle interaction times, these dissipation fields evolve to a relaxed state in 
which they are determined in effect by the dynamical fields and their derivatives. 

A number of technical improvements could be made to strengthen the arguments 
presented here. The physical fluid states for which this result applies are those whose 
gradients are bounded locally to insure that they are not rapidly changing on microscopic 
scales. These local constraints are much stronger than are actually needed to complete the 
proof. All that is really needed are the £ bounds on the fluid fields and their derivatives 
implicit in eqs. (35), (46), and (47). These bounds could undoubtedly be derived using 

far weaker £ 2 conditions on the fluid fields and their derivatives than the local conditions 
used here. A more serious limitation of the present work is its failure to demonstrate 
the existence of any solutions at all of the fluid equations which satisfy these conditions. 
The expectation is that essentially every 'physically relevant' solution to the fluid equations 
does satisfy these conditions. In particular it is expected that 'almost all' initial data which 
are suitably slowly varying on the relevant microscopic length and time scales will evolve 
in such a way that these conditions are preserved for some amount (large on microscopic 
scales) of time. At present, however, theorems of this sort do not exist for these theories. 

Shock waves are one class of physical phenomena that do violate the conditions 
imposed on the fluid states in this work. Significant differences probably do exist in the 
descriptions of this type of fluid phenomenon among the various causal theories and the 
non-causal Navier-Stokes equations. Can meaningful experimental differentiation among 
the various theories be found by observing shock waves? Or, do the predictions of all 
macroscopic fluid theories become meaningless when applied to shocks, since these fluid 
states all contain rapid variations on microscopic particle interaction scales? 
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